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Abstract 



We address a two-dimensional nonlinear elliptic problem with a /mite-amplitude periodic poten- 
tial. For a class of separable symmetric potentials, we study the bifurcation of the first band gap 
in the spectrum of the linear Schrodinger operator and the relevant coupled-mode equations to de- 
scribe this bifurcation. The coupled-mode equations are derived by the rigorous analysis based on 
the Fourier-Bloch decomposition and the Implicit Function Theorem in the space of bounded contin- 
uous functions vanishing at infinity. Persistence of reversible localized solutions, called gap solitons, 
beyond the coupled-mode equations is proved under a non-degeneracy assumption on the kernel of the 
linearization operator. Various branches of reversible localized solutions are classified numerically in 
the framework of the coupled-mode equations and convergence of the approximation error is verified. 
Error estimates on the time-dependent solutions of the Gross-Pitaevskii equation approximated by 
solutions of the coupled-mode equations are obtained for a finite-time interval. 

1 Introduction 

Interplay between nonlinearity and periodicity is the focus of recent studies in different branches of nonlin- 
ear physics and applied mathematics. Physical applications of nonlinear systems with periodic potentials 
range from nonlinear optics, in the dynamics of guided waves in inhomogeneous optical structures and 
photonic crystal lattices, to atomic physics, in the dynamics of Bose-Einstein condensate droplets in 
periodic potentials, and from condensed matter, in Josephson-j unction ladders, to biophysics, in vari- 
ous models of the DNA double strand. The paramount significance for these models is the possibility 
of spatial localization, that is emergence of nonlinear localized structures residing in the spectral band 
gaps of the periodic potentials. To describe this phenomenon, the primary equations of physics are typi- 
cally simplified to the Gross-Pitaevskii equation, which we shall study in our article in the space of two 
dimensions. More precisely, we consider the two-dimensional Gross-Pitaevskii equation in the form 
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where E(x, i) : R 2 xR wC, V 2 = d 2 Xl + d% 2 , V(x) : R 2 ^ R, and a = ±1. Time-periodic solutions of the 
Gross-Pitaevskii equation are found from the solutions of the nonlinear elliptic problem 

V 2 <j>{x) + u<f>(x) = V{x)<p + a\(j)(x)\ 2 <f)(x), (1.2) 

where 4>(x) : R 2 i— > C and w G R arise in the substitution E(x,t) = 4>{x)e~ lujt . It is known that localized 
solutions of the elliptic problem (|1.2j) with a periodic potential V(x), called gap solitons, exist in every 
finite gap of the spectrum of the Schrodinger operator L = — V 2 + V(x) and in the semi-infinite gap for 
a = — 1 [23 US]- Bifurcations of localized solutions from edges of the spectral bands were studied earlier 
inHaEl. 

Coupled-mode equations were used by physicists for the analysis of existence, stability and dynamics 
of gap solitons [18} [26] . A justification of the one-dimensional coupled- mode equations in the context 
of the elliptic problem (|1,2|) with iGM was carried out in [21] in the limit of small- amplitude periodic 
potentials V(x). (A justification of time-dependent coupled-mode equations on finite-time intervals was 
done earlier in [13} 124].) In the limit of small-amplitude potentials, narrow gaps of the spectrum of the 
Schrodinger operator L = —d 2 + V{x) bifurcate from resonant points of the spectrum of Lq = —d 2 , 
while the Bloch modes of L bifurcate from the Fourier modes of Lq. In other words, photonic band 
gaps open generally for small-amplitude one-dimensional periodic potentials V(x). Small- amplitude gap 
solitons of the elliptic problem (|1.2I) in x G R reside in the narrow gaps of L for V(x) ^ according to 
the approximation obtained from the coupled- mode equations |21j . 

We refer to the opening of a spectral gap under a small change of the potential V(x) as to the bifurcation 
of the band gap. Bifurcations of band gaps do not occur for small-amplitude multi-dimensional periodic 
potentials. This is caused by an overlap of the spectral bands of Lq = — V 2 in the first Brillouin zone if 
x G R^ and N > 2 [16]. As a result, photonic band gaps in multi-dimensional potentials open only at 
some finite amplitudes of the periodic potential V{x) and the resonant eigenfunctions are given by the 
Bloch modes of L = —V 2 + V(x) rather than by the Fourier modes of Lq. Although the coupled-mode 
equations were also derived for multi-dimensional problems with small periodic potentials [H El El HI [8] and 
the resonant Fourier modes were used for the approximation of the full solution, the applicability of these 
coupled-mode equations remains an open issue for a rigorous analysis. Bloch mode decomposition has 
been also used in one dimension for finite-amplitude periodic potentials to derive coupled-mode equations 
[9]. The corresponding unperturbed one-dimensional potential, however, has to be of a special type to 
admit a finite number of open gaps, such that a new gap is opened under a small perturbation. 

In this paper we derive coupled-mode equations for wavepackets in narrow band gaps of a finite- 
amplitude periodic potential by using the Fourier-Bloch decomposition and the rigorous analysis based 
on the Implicit Function Theorem in the space of bounded continuous functions vanishing at infinity. The 
coupled-mode equations we derive here take the form of coupled nonlinear Schrodinger (NLS) equations. 
These equations differ from the first-order coupled-mode equations exploited earlier [26] . Similar coupled 
NLS equations have been recently derived in [22] near band edges of the well-separated spectral bands 
and in [6] for tunnelling problems. Unlike these works relying on numerical approximations, we justify 
the derivation of the coupled-mode equations and prove the persistence of localized solutions in the full 
nonlinear problem (II. 2p . Although details of our analysis are given only for the bifurcation of the first 
band gap in the spectrum of L, a similar analysis can be developed for bifurcations of other band gaps 
and for bifurcations of the localized solutions near band edges of the well-separated spectral bands. 

Our derivation is developed for the class of separable potentials 

V(x 1 ,x 2 ) = rj[W(x 1 ) + W(x 2 )], ??GR, (1.3) 
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where the function W(x) is assumed to be real-valued, bounded, piecewise-continuous, and 27r-periodic 
on x £ R. To simplify the details of our analysis, we assume that W(— x) = W{x) on x G R, such that the 
solution set of the elliptic problem (|1.2|) includes functions satisfying one of the following two reversibility 
constraints 



where s\, S2 = dbl. 

Our strategy is to show that there may exist a value rj = r/o, for which the first band gap opens due to 
the resonance of three lowest-order Bloch modes for three spectral bands of the operator L = — S/ 2 + V(x). 
The new small parameter e := rj — 770 is then used for the bifurcation theory of Bloch modes which results 
in the algebraic coupled-mode equations for nonlinear interaction of the three resonant modes. The 
Fourier-Bloch decomposition is used for the approximation of localized solutions and for the derivation 
of the differential coupled-mode equations with the second-order derivative terms. The main idea behind 
our technique is that a differential operator after the Fourier-Bloch decomposition becomes a pseudo- 
differential operator whose symbol is the multiple-valued dispersion relation. Then only the relevant 
(three) branches of the resonant modes can be taken into account, which leads to a coupled-mode system. 

If the linearization operator of the coupled NLS equations at the localized solutions is non-degenerate, 
the persistence of the localized symmetric solutions in the full nonlinear problem (jl.2p is proved with the 
Implicit Function Theorem. Localized symmetric solutions of the coupled NLS equations are approxi- 
mated numerically and the convergence rate for the error of approximation is studied. Finally, we study 
time-dependent localized solutions of the Gross-Pitaevskii equation and the coupled-mode system and 
control smallness of the distance between the two solutions on a finite-time interval. 

The article is structured as follows. Section [2] reviews elements of the Sturm-Liouville theory for the 
separable potentials. Section [3] contains a derivation of the algebraic coupled-mode equations for three 
resonant Bloch modes. Section 2] gives details of the projection technique for the derivation of the dif- 
ferential coupled-mode equations for localized solutions. The persistence of localized reversible solutions 
under a non-degeneracy assumption on the linearization operator is proved in Section [5j Numerical ap- 
proximations of the localized solutions, the associated linearization operators, and the convergence of the 
approximation error are obtained in Section [U Section [7] extends the results to the time-dependent case 
for finite time intervals. Section [8] discusses relevant generalizations. 

2 Sturm— Liouville theory for separable potentials 

It is typically expected that the band gaps in the spectrum of the linear Schrddinger operator with a 
two-dimensional periodic potential open at the extremal values of the Bloch band surfaces |16j . For a 
general potential, however, the extremal values may occur anywhere within the first irreducible Brillouin 
zone Bq in the quasi-momentum space (ki,k 2 ) [Hj. We will show here that the extremal values for a 
separable potential occur only at the vertex points T, X and M on the boundary 8Bq. Figure Q] shows 
the irreducible Brillouin zone Bq and the vertex points for a two-dimensional separable potential V(x). 

Let us consider the spectral problem for the linear Schrddinger operator associated with the separable 
periodic potential f)l .3|) : 




(1.4) 



or 




(1.5) 



V 2 u(xi,x 2 ) + rj [W(xi) + W(x 2 )] u(xi,x 2 ) = oju{xi,x 2 ). 
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Figure 1: The first irreducible Brillouin zone Bq for a two-dimensional separable potential. 

By using the separation of variables u(x\,X2) = fi{xi)f2(x2) and the parametrization to = £1 + we 
obtain the uncoupled eigenvalue problems 

- f!{xj) + n \Y {.,■,))■ ,{.,■,) = Cjfjixj), j = 1, 2. (2.2) 

The Bloch modes and the band surfaces for the one-dimensional problems (|2.2j) are introduced according 
to the regular Sturm-Liouville problem 



—u"{x) + r]W(x)u(x) = pu(x), < x < 2ir, 
u(2tt) = e i27rfc u(0), 



/"2-A-,,/n\ (2-3) 



where A: is a quasi-momentum defined on the interval T = [— ^, ^] • The eigenfunctions of the Sturm- 
Liouville problem (|2.3j) are periodic with respect to k with the period one. By Theorem 2.4.3 in |10j . 
there exists a countable infinite set of eigenvalues {p n (k)} n en for each k € T, which can be ordered as 

Pi(k) < P2 (k) < P3 (k) < ... 

By Theorem 4.2.3 in [ID], if W(x) is a bounded and 27r-periodic potential, eigenvalues {p n (fc)}neN have 
a uniform asymptotic distribution on k £ T, such that 

C_n 2 < \p n {k)\ < C+n 2 , VneN, Vfc £ T, (2.4) 

for some constants C± > 0. 



Let {u n (x;k)}[£^ be the corresponding set of eigenfunctions of the Sturm-Liouville problem (|2.3p . 
such that n n (a; + 2n;k) = u n (x; k)e 27rik . By Theorem XIII.89 in [23], if W(x) is a bounded, piecewise- 
continuous and 2-7r-periodic potential, then the eigenvalue p n (k) and the Bloch function u n (x;k) are 
analytic in k G T\ {0, and continuous at the points k = and = By Theorem XIII. 95 in [23], 
the eigenvalue p n {k) is extended on a smooth Riemann surface in the neighborhood of the points k = 
and k = ^ if the n th spectral band is disjoint from the adjacent (n ± l) th spectral bands by a non-empty 
band gap. 

By Theorem XIII. 90 in |23j, if W{x) is a bounded, piecewise-continuous and 27r-periodic potential, 
then the spectrum of Lid = —d 2 + rjW(x) in L 2 (M) is absolutely continuous and consists of the union of 
the intervals in the range of the functions p n {k) on k £ T for n £ N, where p n {—k) = p n {k). Moreover, 
the extremal values of p n (k) occur only at the points k = and k = ±k. These points correspond to an 
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alternating sequence of the maximum and minimum values of the eigenvalues {p n (k)} n £N according to 
the following formula 

argmin/)2m-i(^) = 0, argmaxp 2 m-iO) = ±7:, aig min p 2 m (&) = ±7:, argmax p 2m (k) = 0, 
fcss Kg§ z fce§ z fees 

for all m G N. The Bloch function u n (x; k) is 27r-periodic if k = and 27r-antiperiodic if k = ±5. Let us 
denote the periodic and antiperiodic eigenfunctions and the corresponding eigenvalues by 

i> n (x) = u n (x; 0), A n = p n (0) and (p n (x) = u n (^x; , p n = p n (±7^) > n G N - ( 2 - 5 ) 

By Theorems 2.3.1 and 3.1.2 in [10], these eigenvalues are ordered by Ai < A2 < A3 < ... and p± < P2 < 
P3 < •••) while the corresponding eigenfunctions ip n {x) and <p n (x) have precisely n — 1 zeros (nodes) on 
the interval (— 7r, 7t). If W(— x) = W(x), the eigenfunctions tj} n {x) and tp n {x) are even for odd n and odd 
for even n. Each set of the eigenfunctions {tpn{ x )}neZ and {<p n (z)}n&z is orthogonal in L 2 ([— ir, tt]). 

Eigenvalues of the same one-dimensional operator L\d = —d% + rjW(x) for 4-7r-periodic eigenfunctions 
consist of the union of the eigenvalues {A n } n6 N and {// n } ng N- Since there are at most two eigenvalues of 
the second-order operator L\d and by Theorem 2.3.1 of [10], we obtain that 



Ai < pi < P2 < A2 < A 3 < /x 3 < pa < A 4 < A 5 < ... (2.6) 

In particular, the first band gap is always non-empty for a non-constant potential W(x) (see Theorem 
XIII. 91(a) in |23j). Due to this ordering, the lowest value of u = p ni (k\) + p n2 (hz) for the crossing 
of the Bloch band surfaces associated with the two-dimensional separable potential V(x\,X2) occurs at 
uj = u>o = Ai + P2 = 2pi- Using these facts, we obtain the following results. 

Lemma 1 Extremal values of the Bloch band surfaces to = p ni (k{) + p n2 {k2) for (ki, k^j G Bq occur only 
at the vertex points of the boundary 8Bq . 

Proof. Because uj = p ni {k\) + p n2 (^2) for a separable potential (|1.3I) . we have Vw = [p' ni (ki), p' n2 (k2)] T , 
where p n (k) are analytic on /c G T\ |0, ±^}. If an extremal point occurs in the interior of the Brillouin 
zone Bq, then Vw = at k = ko G -Bo- However, p' n {k) / for any < |fc| < ^ and any n G N by 
Theorem XIII. 90 in |23| . Similarly, the extremal points cannot occur in the interior of the boundary 8Bq. 
Therefore, the extremal values of uj occur only at the vertex points T, X and M on 8Bq. □ 

Remark 1 The derivative p' n (k) may be non-zero at k = and k = i if the n th spectral band touches 
the adjacent (n ± l) th spectral bands. This happens when the corresponding eigenvalue A n or p n is 
double degenerate with the equality sign in the ordering (|2,6p . However, Ai, pi and P2 are always 
simple and, therefore, Vuj = at least for the first three spectral bands uj = p ni (ki) + p^ifo) at 
uj = uj = Ai + P2 = 2pi. 

Lemma 2 Assume that the resonant condition \\ + P2 = 2p\ = ujq for the 2ir-periodic eigenvalue Ai 
and the 2tt- antiperiodic eigenvalues p\,p2 of Lid = — d% + r]W(x) is satisfied for i] = t]q. Then there are 
exactly three resonant Bloch modes at uj = ujq and r] = ijq in the spectral problem \2. 

$1 = "01 (^1)^2(^2), $2 = ( P2{x\)i>l(x2), $3 = (pi(xi)ipi(x 2 )- 
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The three resonant modes are orthogonal to each other with respect to the Air-periodic inner product 

/2tv /*2-7r 
/ f(x 1 ,x 2 )g(xi,x 2 )dx 1 dx 2 . 
-2tt J-2-k 

Proof. There are exactly three resonant modes for any separable potential (|1,3|) because Ai and [i\ are 
the smallest non-degenerate eigenvalues of the operator L\d and the double degeneracy Ai + fi 2 = + Ai 
is due to the symmetry with respect to the interchange of the variables x\ and x 2 . The eigenfunc- 
tions {ip n (x)} n ^ and {tp n {x)}neN are an orthogonal to each other in L 2 ([— 2tt, 2tt]), which leads to the 
orthogonality of the three resonant modes. □ 

Remark 2 We note that the eigenfunctions ipi(x) and tpi(x) are not orthogonal on the interval [— vr,7r] 
since both of them are positive by Theorem 3.1.2 in [10J. However, a linear combination of ipi(x) and 
(pi(x) belongs to the class of 47r-periodic functions and the two eigenfunctions are orthogonal on the 
double-length interval [— 2ir, 2tt]. 

Remark 3 By Theorem 6.10.5 in [TO], there are finitely many gaps in the spectral problem (|2.ip with 
a separable potential. (In fact, there are finitely many gaps for general smooth periodic potentials in 
dimension 2 or higher [20J.) Lemma [2] indicates a bifurcation when one of these gaps at the lowest value 
of u) may open. It is observed for many examples of the separable potential (jl.3p that this bifurcation 
leads to the first band gap in the spectrum of L = — V 2 + 77[W(xi) + W(x2)] if f] is increased from r\ = 0. 

Example 1 Let W(x) = 1 — cosx. Numerical approximations of the eigenvalues of the Sturm-Liouville 
problem (|2.3p are computed with the use of the second-order central difference method. The eigenfunc- 
tions ip n (x) and <p n (x) are plotted in Figure[2]for n = 1,2, 3, while the dependence of the first eigenvalues 
A n and fj, n on rj is plotted in Figure EJa). Figure [3jb) shows the intersection of 2/ii and Ai + fi 2 at 
the bifurcation value r/ = rjo « 0.1745, when Ai « 0.1595, fj,% rj 0.3336 and fi 2 0.5077, such that 
<^o = Ai + Hi = 2/ii 0.6672. The perturbation behavior shown in Figure [3] was studied analytically in 
[5j. Figure [5] illustrates the band structure along OBq of the full spectral problem (|2.ip with r/ = 0.1745 
clearly revealing the resonance at oj ~ 0.6672 and the two resonant modes at the points X and M. The 
third resonant mode lies at X' because the symmetry V{x\, x 2 ) = V(x 2 , x\) implies u)(ki, k 2 ) = ui(k 2 , k\). 

In what follows, we consider the bifurcation of non-trivial solutions of the nonlinear elliptic problem 
(jl.2p with the separable potential (jl.3p in the lowest band gap described by Lemma [2l Let e = r/ — rjQ, 
u) = ujq + e$7, <p(x) = y/e<&(x) and rewrite the nonlinear elliptic problem (jl.2p in the form 

C <S>{x) = efl^{x) - e [W(x x ) + W(x 2 )] - eo-\<S>(x)\ 2 $(x), (2.7) 

where $ : M 2 i— ► C and Cq = —V 2 + tjq [W(x\) + W(x2)] — ujq. Two classes of non-trivial solutions of 
the bifurcation problem (|2.7p are considered for small e: bounded 4-7r-periodic solutions (Section [3]) and 
bounded decaying solutions (Sections HHS]). 

3 Algebraic coupled-mode equations 

We consider here bounded periodic solutions <3?(x) of the bifurcation problem (|2.7p for small e. Our 
results depend on the period of the periodic solutions $(x). Since the potential V(x) is separable, the 
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Figure 3: (a) Eigenvalues of Lid = — d^. + r](l — cosx) for 27r-periodic (solid) and 27r-antiperiodic (dashed) 
boundary conditions, (b) The lowest resonance occurs at lu = ujq and r\ = r/o when \\ + ^2 = 2//i. 
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Figure 4: Band diagram along 8Bq for the spectral problem (|2.ip with W(x) = 1 — cosx and r] = t]q. 

27r-periodic (27r-antiperiodic) function <&(x) can be represented by the series of Bloch modes associated 
with the eigenvalue problem f|2.3[) on x G R for k = (k = ±^), while the 4-7r-periodic functions 
can be represented by the series of both 2-7r-periodic and 27r-antiperiodic Bloch modes. In what follows, 
we shall use symbols L 2 (Pfc), H s (Fk) and CP(Pfc) to denote the corresponding spaces of functions u(x) 
on the interval [0, 27r] satisfying the boundary conditions u(2ir) = e l27Tk u(0). We shall also use symbol 
to denote the vector of elements of the sequence {4> n } n< =^. 

Proposition 1 Let W(x) be a bounded and 2-n-periodic function. Let {p n (k)} n eN and {u n (x; A;)} ng N be 
sets of eigenvalues and eigenf unctions of the Sturm-Liouville problem h2.3\) that depends on k E T, such 
that 

{u n (-;k),u n >(-;k)) 2 TT = 5n,n', Vn,n'eN, \/k G T, (3.1) 

where (/, 5)2^ = f{x)g{x)dx and 5 n , n i is the Kronecker symbol. For any fixed k G T, the set of 
eigenf unctions is complete in L 2 (Pfc) ; such that there exists a unique set of coefficients {4> n }neN in the 
decomposition 

n u n (x, k), (3.2) 

neN 

given by <f) n = (u n (-;k),(f))2 7 r- 

Proof. The statement of the proposition follows by Theorem XIII. 88 of [23j . □ 

Lemma 3 Let 4>(x) be defined by the decomposition Ii3.^) . where d> belongs to the vector space l](N) with 
the norm 

n\k(N) = Y,( i+n y^ <o °- ( 3 - 3 ) 

ngN 

If s > 2> then 4>(x) is a continuous function oni£l and 4>(x) = r ip(x)e lkx with tp(x + 2-zr) = ip(x). 
Proof. By the triangle inequality, we obtain 

IHIc°(R) < ^2 I^IIKO;&)llc°(p & )- 
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By Sobolev's Embedding Theorem, there exists a C > such that ||u„(-; &)llc (P*.) — ^\\ u n('i k)\\H s (¥ k ) 
for any s > |, n G N and i£l Since VK(x) is a bounded potential, the squared norm ||ii n (-; &)||#s(p fe ) 
is equivalent to the integral 

(c„ + L 1D ) s/2 u n (x; k) dx = (c„ + p n (k)) s , 

where c v > — r/ min^gR W(x) and (c n + L\d) s I 2 is defined using the spectral family associated with the 
Sturm-Liouville operator Lid = —d 2 + r]W(x) in L 2 (Pfc). By the asymptotic distribution of eigenvalues 
(|2.4p . we obtain that 



OS 



\) ^ C Y1 \<M\M';k)\\H'(p k ) <c^{i + n) s \<i) n \ 



b 1 

nSN neN 

for some C > and any s > ^. □ 

If the potential r]W{x) is continued with respect to the parameter rj, the perturbation theory for 
an eigenvalue p n (k) depends on whether p n {k) is simple or multiple. The following lemma is a trivial 
statement of the perturbation theory, so we omit its proof. 

Lemma 4 Let W{x) be a bounded and In -periodic function of x Gl. Fix k £ T ; r/o € R, and n G N. // 
Pn-i{k) < p n (k) < p n+ \{k), then p n {k) and u n (x;k) depend analytically on rj near r] = rjo, such that 

d v Pn(k)\ v = V o = {un(-;k)\ v=V0 ,Wu n (-;k)\ ri=V0 )2n (3.4) 

and 

\Pn(k) - p n {k)\ v=V0 - d v p n (k)\ v=710 {ri - r) )\ < C(r] - i] ) 2 (3.5) 
for some C > and sufficiently small \r] — rjo\- 



Due to the construction of the spectrum of the two-dimensional spectral problem (|2.ip . the results of 
Lemma H] settle the question of the splitting of the three resonant bands near the point to = ujq described 
in Lemma [2j Indeed, each Bloch band surface (the graph of u> = p ni (k\) + p n2 (/C2)) corresponding to the 
three resonant bands includes a different vertex point {X' , X, and M) at u = ujq for different values of 
k in the set {(0, ^) ; (^,0) ; Q, ^)}. Each eigenvalue of the spectral problem (|2.ip is simple and isolated 
if k is fixed, such that the continuation of each Bloch band surface in 77 near r} = 7/0 is determined by 
the expression (13. 4p and the separation of variables in uj = p ni {k\) + pn^ik^). In other words, bifurcation 
of the lowest band gap that occurs at uj = ujq can not be studied if the function space is restricted to 
L 2 (Pfc 1 x Pfc 2 ) for a fixed value of (k±, k2) £ T 2 . To detect the bifurcation, we shall work in the function 
space L 2 ei ([— 2ir, 2ir] x [— 27T, 2-k]), where all three vertex points X' , X and M correspond to the same 
47r-periodic boundary conditions on the function <I>(x). According to Lemma [21 the three resonant Bloch 
modes are orthogonal to each other in this space. 

To work with the 4-7r-periodic eigenfunctions, we shall introduce a number of notations. Let {^ n }ngN 
and {u n (x)} n£ N denote the sets of eigenvalues and the corresponding 47r-periodic eigenfunctions of the 
one-dimensional spectral problem (|2.3p . These sets consist of the union of eigenvalues and eigenfunctions 
defined by (|2,5p . Suppose that the eigenvalues are sorted in non-decreasing order, such that v\ = Ai, 
v 2 = pi-, V3 = P2-, ^4 = A2, = A3 and so on, and that the eigenfunctions v\ = tf)\,V2 = <pi,vs = 
<P2,V4 = ^2,^5 = etc., are normalized such that (^,1^)4^ = 5ij, where {f,g)A_-K = J(x)g(x)dx. 
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The statements of Proposition [T] and Lemma [3] extend naturally to the new sets of eigenfunctions in space 
Lp er ([— 2-7T, 27r]). To develop the nonlinear analysis of bifurcations of a non-trivial 47r-periodic solution 
<&(x) on i 6 R 2 , we shall study first the nonlinear vector field |<&(a;)| 2 $(a;) acting on the decomposition 
&(x) = SneN 4>nV n {x)i where <j> belongs to the vector space iJ(N). 

Lemma 5 LetW(x) be a bounded and 2ir -periodic function. Let(j)(x) = XlneN ( t > n v n{ x ) an d |<^>(a;)| 2 <^(a;) = 
J2neN 9n v n{x) . If < s < 1, then there exists a C > such that ||g||ji(N) — ^[|0llfi(N) • 

Proof. We only need to prove that the space /](N) with < s < 1 forms a Banach algebra in the sense 

V4>,tp G Z*(N) : ||0*0lbi(N) < C|l^ll«i(N)ll^ll«i(N), (3-6) 

for some C > 0, where 

(4>*<p) =^2^2 K n,i,j(f>Wj, VneN 

and 

#n,i,j = (v n ,ViVj)4n, V(n,i,j) e N 3 . 
According to Lemmas A.l and A. 2 of Appendix A of [7J, if W € Lp er ([— 27r, 2tt]) and ||v n ||f^ ^ 2lT 2irl) = ^ 
for any n S N, then there exists an n-independent constant C > 0, such that 

< tt— I C . rrTT, V(n,i, j) e N 3 , (3.7) 

(1 + |n — z — j | J p 

for any < p < 2. By explicit computation, we obtain 

\\4>*<P\\ll(N) < 53(l + ») < ' J^I-K'n.ijIl^illVil 

< cya+iri^iya+jri^iy f 7 — ^ — r V 7 — — 1 — iT - 

ieN jeN neN vv A J > J K 1 J|; 

< cy^n+iYUA y(i+j) s iw 7 -i y f i + - — -^7 — , 1 r 

j£N jGN neN V V ' J \ J J / \ I 

for some C > 0. If jj > 1 and p — s > 1, the bound is completed as follows 

H0*0l|ji(N) < C'l|l0ll/](N)ll ( ?ll/](N) + ^2||0||ji(N)||0||ii(N) < {pi + C2) ||^||ji(N)ll0llii(N) 

for some Ci,C 2 > 0. Since 1 < p < 2, the parameter s must satisfy < s < 1. □ 
By using this construction, we prove the first main result of our analysis. 

Theorem 1 Let W(x) be a bounded, even and 2ir-periodic function on x £ R. T/ie nonlinear elliptic 
problem \2. 7\ ) has a continuous and At: -periodic solution <&(x) for sufficiently small e if there exists a 
solution for (A\, A 2: A 3 ) G C 3 of the algebraic coupled-mode equations 



(Q-/?iL4i = a 
(Q-Pi)A 2 = a 
(n-f3 2 )A 3 = a 



11 \A 1 \ 2 A 1 +- f2 {2\A 2 \ 2 A 1 +A 2 A 1 )+- f3 (2\A 3 \ 2 A 1 +A 2 A 1 ) 
-f 1 \A 2 \ 2 A 2 + -f 2 (2\A 1 \ 2 A 2 + AjA 2 ) +-/ 3 (2\A 3 \ 2 A 2 + AjA 2 ) 



+ eR 1 (A 1 ,A 2 ,A 3 ), 
+ eR 2 (A 1 ,A 2 ,A 3 ), 



74 |^ 3 | 2 ^3 + 2 73 (|Ai| 2 + \A 2 \ 2 )A 3 + l3 (A 2 + A 2 )A 3 ] + eR 3 (A 1 ,A 2 ,A 3 ), 

(3.8) 
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where parameters are given by 

0X = (l/>l,Wlpx)4K + (<P2,Wip 2 )4Tr, #2 = 2{if U W^l)^ 

and 

71 = (^1,^1)477(^2^2)4^, 72 = K^i, ^1>4tt| 2 , 73 = (^1, ^1)4^1, ¥>i)47T, 74 = |(^l,^l)47r| 2 

and the residual terms ii^^^^i, ^2, ^3) are analytic functions of e near e = satisfying the bounds 

V|e| < e : \R 1) 2,3(A 1 ,A 2 ,A 3 )\ < Ci,2, 3 (|^i| + \A 2 \ + \A 3 \), 

for some constants C123 > and sufficiently small €q. Moreover, there exists an e-independent constant 
C > such that 

||$ - Ax^x - A 2 $ 2 - A 3 <S> 3 \\ c o {R2) < Ce, 
where (<l?i, $2, $3) are the modes defined by Lemma\M 

Proof. 4ir— periodic solutions of the nonlinear elliptic problem (|2.7p are expanded in the form 

${x)= ^ ^n 1 ,n 2 V ni (xi)v n2 (x 2 ). (3.9) 
(ni,n 2 )eN 2 

Let the vector <1> with the elements of the sequence {$ni,n 2 }(ni.n 2 )6N 2 belong to the vector space Z](N 2 ) 
equipped with the norm 

ll*lbl(N 2 )= Yl (l + ^l) S ( 1 + ^2)1^n 1 ,n 2 | < OO. (3.10) 

(ni,n 2 )eN 2 

Similarly to Lemma [3l it follows that the function $(x) is continuous if $ G ^(^ 2 ) with s > 5. Indeed, 

ll^llc 6 °(M 2 ) - I < ^ni i «2lll' u nillc0([-27r,27r])ll t 'n 2 llc0([-27r,27r]) 

(m,n 2 )eN 2 

< C |^ > n 1 ,n 2 |||'V„ 1 ||Hs cr ([-27r,27r])ll t 'n 2 ll^ cr ([-27r,27r]) 

(m,n 2 )eN 2 

< C (l + ni) s (l + n 2 ) s \$ ni , n2 \ 

("i,n 2 )eN 2 

for some C, C > and any s > ^. By construction, the solution <£(x) is 47r-periodic in both coordinates 
of x G R 2 . The partial differential equation (|2.7p is rewritten in the lattice form, which is diagonal with 
respect to the linear terms 

Kj + v n2 ~ uj ~ efl] = -ea ^ M n , m4J <!> m $>i$> j , Vn = (ni,ra 2 ) G N 2 , (3.11) 

(m,ij')eN 6 

where v n depend on r\ = tjq + e and 

,m,i,j ("f ni Vii , "mi ^ji )47r ("^712^12 > ^12^,72) 47T, V(n 1 ,n 2 ,m 1 ,m 2 ,ii,i 2 ,ji,j 2 ) G N 8 . 
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Since M n>m ij is a product of one-dimensional inner products and the weights in the norm in ^(N 2 ) are 
separable, Lemma [5] applies and guarantees that the nonlinear vector field of the lattice equations (|3.1ip 
is closed in Z](N 2 ) for < s < 1. Therefore, solutions of the lattice equations (|3.1ip can be considered in 
the space ^(f^ 2 ) f° r an y | < s < 1. 

By Lemma [H the three resonant modes are isolated from all other 4-7r-periodic modes. By Lemma [H 
the eigenvalues v ni + v n2 are then analytic in rj near rj = tjq. The three resonant modes correspond to 
the values of n in the set {(1, 3); (3, 1); (2, 2)}. Therefore, we decompose the solution $ of the lattice 
equations (|3.1ip into 

$ = Axe^ + A 2 e 3A + A 3 e 2 ,2 + & , (3.12) 
where {ex,3, 63,1, 62,2} are unit vectors on N 2 , Span(ei j 3, e^ t i, 62,2) is the kernel of the linearized system 
at the zero solution for e = 0, and \& lies in the orthogonal complement of the kernel such that ^1,3 = 
^3,1 = ^22 = 0. The linearized operator projected onto the orthogonal complement of the kernel is 
continuously invertible for sufficiently small e. By the Implicit Function Theorem in the space ^(N 2 ) for 
any ^ < s < 1, there exists a unique map \& e (Ai, A 2 , A3) : C 3 1— » /g(^ 2 ) f° r sufficiently small e. Moreover, 
the map is locally analytic in e near e = 0, such that ^o(A\, A2, A3) = 0. In addition, ^(0,0,0) = 
for any eel and (d e Sfr e (Ai, A 2 , A 3 ))\ e =o is a homogeneous cubic polynomial of (Ai, A 2 , A3). Let S = 
\A\\ + \A 2 \ + I -^3 1 < ^0 for a fixed e-independent 5o > 0. Then the map & e (Ai, A2, A3) satisfies the bound 

V|e|<e : ||* e (Ai, A 2 , A 3 )\\ lim < eC {\A X \ + \A 2 \ + \A 3 \) , (3.13) 

where eo > is sufficiently small, 5q > is finite, and the constant C > is independent of e and 5. 

We can now consider the three equations of the system (|3.1ip for n = {(1, 3); (3, 1); (2, 2)}. After 
the Taylor expansion of (|3.1ip in e the formula (|3.4p yields the coefficients f3\ , /3 2 since = d e . Note 
that the 27r-inner product in (|3.4p is replaced by the 47r-inner product. Although the Bloch modes 
v n (x) are now normalized over the 4-7r-long interval, the value of (|3.4|) is unchanged. Using the bound 
(|3.13p and the fact that the lattice equations are closed in Z^(N 2 ), we derive the algebraic coupled-mode 
equations in Theorem [TJ The coefficients of these equations can be simplified due to the fact that 
W(x) is even on x G K, such that the eigenfunctions -01 0*0 and <fi(x) are even while <f2(%) is odd on 
x £ R. As a result, many coefficients of the coupled-mode system are identically zero, for instance 
{'4>i l P2, l pi)4iT = 0, (ifi, </?iV?2)4tt = and so on. The residual terms -Ri, 2,3(^1, A2, A3) are estimated from 
the map \& e (Ai, A 2 , A3) with the bound (|3.13p . The last estimate in the statement of Theorem [T] is 
obtained from f)3. 13|) and Lemma [3j □ 

We shall refer to the system (j3.8j) without remainder terms tRi,2,3{Ai, A 2 , A3) as to the truncated 
coupled-mode system. 

Corollary 1 There exist five invariant reductions of the truncated coupled-mode system, namely (i) 
A\ = A2 = 0, (ii) A\ = A3 = 0, (Hi) A2 = A3 = 0, (iv) A\ = 0, and (v) A2 = 0, which persist in 
the full lattice equations \3.11\) . 

Proof. Any of the five invariant reductions implies symmetry constraints on the function <3?(x) at the 
leading order of the decomposition (|3.12p . For instance, if A\ = A2 = 0, then <&{x) is 2-7r-antiperiodic 
with respect to both x\ and x 2 ; if A\ = 0, then <J?(x) is 2-7r-antiperiodic in x\ and 4-7r-periodic in x 2 , and 
so on. By the completeness results of Proposition [U all other terms of the decomposition (13. 9p which 
violate the symmetry constraints on the solution <l?(x) can be set to be identically zero. By the Implicit 
Function Theorem, the zero solution is unique near e = 0. Therefore, the series (|3.9p shrinks to fewer 
terms, the reduction persists for sufficiently small e, and the proof of Theorem [1] applies. □ 
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Remark 4 The invariant reduction A3 = of the truncated coupled-mode system may not satisfy the 
full lattice equations (|3.1ip since the solution with A3 = is still a 47r-periodic function of both x\ 
and X2 and the series (|3.9p can not be shrunk to fewer terms for A% = 0. 

Remark 5 The Fourier-Bloch decomposition needed in Theorem [T] can be alternatively developed for 
in the vector space if(N) equipped with the squared norm 

ii^ii|(N) = E( 1+n2 ) s i^i 2<00 - ( 3 - 14 ) 

neN 

Indeed, the squared norm H^H^flPj,) is equivalent to the integral 

dx = ^2 ^ n i^ n 2 ( c v + Pn(k)) s (u ni (-; k), M„ 2 (■;k)) 27r 

ni€N H2GN 

= ^2 ( C V + Pn(. k )Y l^n] 2 , 
neN 

where we have used the orthogonality relation (|3.ip . By the asymptotic distribution (|2.4p . [|0||^»np fc ) ^ s 
thus equivalent to ||0|| 2 2( N )- By Sobolev's Embedding Theorem, ||</>|lc°(P fc ) — C|l^ll-ff s (P fc ) f° r some C > 
and any s > 5- Therefore, the decomposition (|3.2|) produces a continuous function <j>{x) on x G R if 
4> G ig(N) with s > \. Furthermore, using Lemma 3.4 in [7], one can prove that the nonlinear term maps 
£ 2 (N) to / 2 (N) for any s > 5, such that the arguments of the Lyapunov-Schmidt reductions of Theorem 
[1] will work in the space Z 2 (N) for any s > \ . Note that this approach would lift the upper bound on the 
index s in Theorem [H where \ < s < 1. 

Example 2 For W(x) = 1 — cos a; and r\ = r/o ~ 0.1745, the parameters of the algebraic coupled-mode 
equations (|3.8p are approximated numerically as follows: 

fa « 2.2835, /3 2 « 0.9183 

and 

71 « 9.4829 x 10~ 3 , 72 « 4.5196 x 10~ 3 , 73 « 3.7942 x 10~ 3 , 74 « 1.5981 x 10~ 2 . 

The algebraic coupled-mode equations describe both linear and nonlinear corrections to the eigenvalues 
u = Pm(ki) + ^2(^2) f° r fixed values of A^) at the vertex points X, X' and M. In particular, these 
corrections show how the band edges of the three resonant bands split for e / and deform due to 
nonlinear interactions between resonant Bloch modes. We skip further analysis of the algebraic coupled- 
mode equations and proceed with the decaying solutions of the nonlinear elliptic problem (|2.7p described 
by the differential coupled- mode equations. 

4 Differential coupled-mode equations 

We consider here bounded and decaying solutions <5(x) of the bifurcation problem (12. 7p for small e. The 
decomposition of the decaying solution 3>(x) depends now on the completeness of the Bloch modes of the 
one-dimensional Sturm-Liouville problem (12. 3p in L 2 (M), where all Bloch modes for all k G T must be 
incorporated in the Fourier-Bloch decomposition. 



( C „ + L 1D y/ 2 <p(x) 
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Proposition 2 Let W{x) be a bounded and In -'periodic function. Let {p n (k)}neN an d {u n (x;k)} ne ^ be 
sets of eigenvalues and eigenf unctions of the Sturm-Liouville problem \2.3\) on k £ T, such that 

{u n (-;k),u n >(-;k')) R = 5 njn >5(k - k'), Vn,n'eN, V/s, k' G T, (4.1) 

where {f,g)m = J^f(x)g(x)dx and 5(x) is the Dirac's delta function in the distribution sense. Then, 
there exists a unitary transformation T : L 2 (R) i— ► Z 2 (N, L 2 (T)) given by 

G L 2 (R) : 0(fc) = T<f>, 4> n (k) = [ u n (y; k)^y)dy, Vn G N, V/c G T, (4.2) 
where Z 2 (N, L 2 (T)) is equipped with the squared norm 

IMI/ 2 2(N,L2(T)) = E / \4>n(k)\ 2 dk. 

The inverse transformation is given by 

\/(f> G / 2 (N, L 2 (T)) : <f)(x) = T~ 1 ^ = V [ ^ n (k)u n (x;k)dk, Vx £ R. (4.3) 



Proof. The original proof of the proposition can be found in [11] . Orthogonality and completeness of 
the Bloch wave functions {u n (x; k)} n ^ on k G T is summarized in Theorem XIII. 98 on p. 304 in [23]. 
Existence of the unitary transformation (|4.2p with the inverse transformation (|4.3p is summarized in 
Theorems XIII.97 and XIII. 98 on pp. 303-304 in [23]. □ 

Lemma 6 Let <^>(x) be defined by the decomposition j^.3| ), where <p belongs to the vector space Z](N, L l (T)) 
with the norm 

"'zi(N,Li(l 

neN neN 



(T)) = E^ 1 + n YUn\\v(j) = £(1 + n ) S / \k(k)\dk < oo. (4.4) 

neN neN ^ T 

If s > h, then (f)(x) is a continuous function on x G K sitc/i i/iai ^>(x) -» as |x| — * oo. 



Proof. The proof is similar to that of Lemma [3] since the asymptotic bound (|2.4p is uniform on fc G T. 
Therefore, 



< 

neN " R neN 



E / \^ri{k)\\\u n {-,k)\\ c o {¥k) dk < / l0n( fc )llkn(-; fc )ll^(P fe )^ < C'll'Alli](N,Li(T)) 



for some C, C > and any s > ^. The decay of 4>(x) as |x| — > oo follows from the Riemann-Lebesgue 
Lemma applied to the Fourier-Bloch transform. Indeed, the summation in n G N of the integrals on 
k G T can be written as an integral on k G M. in the form 

4>(x) = / (j)(k)v(x; k)dk, 
Jr 

where v(x;k) for k G M is related to the Bloch functions u n (x;k) for k £ T (see [7] for details). Since 
v(x;k) = e lkx w(x; k), where w(x;k) is periodic in x with period 2tt and is uniformly bounded in C^(1R) 
with respect to both x and k, the Riemann-Lebesque Lemma applies to the Fourier-Bloch transform in 
the same manner as it applies to the standard Fourier transform if G L 1 (M). □ 
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Lemma 7 Let W{x) be a bounded, piecewise- continuous and 2ir -periodic function. Let p n {ko) be the 
extremal value of p n (k) for either ko = or k$ = ±5 and assume that the adjacent bands are bounded 
away from the value p n {ko). Then, p n (k) is an analytic function at the point k = k$, such that 



P '(k ) = o, pM^o, 

for some C > uniformly in k G T. 



Pn(k) - p n (k ) - ^Pn(k )(k - k ) 2 



< C(k - k ) 4 (4.5) 



Proof. By Theorem XIII. 95 in [23j, the function p n (k) for an isolated spectral band can be extended onto 
a smooth Riemann surface which has no singularities other than square root branch points at k = and 
k = ±i. Therefore, the function p n (k) is expanded in even powers of (k — ko) for ko = or ko = ±^ and 

P'n(k ) + 0. □ 

Lemma 8 Let W(x) be a bounded and 2tt -periodic function. Let (j) = Tcj) and g = T (\(j)\ 2 <p) . IfO < s < 1, 
then there exists a C > such that WgWiKfn^ 1 ^)) < C||0llfi(N l^t)) ■ 

Proof. The proof is similar to that of Lemma [5] since the asymptotic bound (|3.7p is uniform on k G T 
|7j and the length of T is finite. The Banach algebra property (13.61) holds in space /](N,L 1 (T)) for any 
< s < 1. Indeed, let 



[4>*<p) (k) = V] y2 / K ndl j 2 (k,ki,k 2 )(f>j 1 (ki)(p j2 (k2)dkidk2, n G N, fe G T, 

where if njl j 2 (A;, fei, k 2 ) = (%(•; fe), «ji (•; h)u j2 (-; k 2 )) R i for all (n,jx,j 2 ) G N 3 and (k,ki,k 2 ) G T 3 . By 
an explicit computation, we obtain 

H#*0lllJCN,£i(T)) < E^ + ^E E / / \ K nj u j 2 (k,k 1 ,k2)\\<p jl (k 1 )\\ip j2 (k 2 )\dkdk 1 dk2 

jieN j 2 eN 



'T JT 

X ^V(i + ii)(i + i 2 )y (l + |n-. 



neN 



31 ~J2 



for some C > and any < p < 2. The last sum is finite if p > 1 and p — s > 1 (that is < s < 1), 
similarly to the proof of Lemma □ 

We shall now define the coupled-mode system, which is a central element of our paper. Let Ai 2 3(2/) 
be functions of y = y^x on 1 6 I 2 which satisfy the differential coupled-mode equations 



y 1 \A 1 \ 2 A 1 + 72 (2|A 2 | 2 A 1 + Ap x ) + 73 (2|A3| 2 A + Ap! 
7i|^ 2 | 2 ^2 + 7 2 (2|Ai| 2 A 2 + AjA 2 ) + 7 3 (2|^ 3 | 2 ^ 2 + A\A 2 
74 |^ 3 | 2 ^3 + 273d A| 2 + \A 2 \ 2 )A 3 + 73(^4? + A 2 )A 3 ] , 

(4.6) 

where parameters /3i 2 and 71,2,3,4 are the same as in Theorem [IJ while 



(n-(3 1 )A 1 + (a 1 d 2 yi +a 2 d 2 jA 1 = a 
(n-P 1 )A 2 + (a 2 d 2 1 +a 1 d 2 2 )A 2 = a 
(n-p 2 )A 3 + a 3 {d 2 1 +d 2 2 )A 3 = a 



' 1 I = ^Pl(O)) "2 = ^P 2 Q^) , " -//, | - 



1 « 



Since Ai = pi(0), /«i = pi (|), and p 2 = P2 (3) are non-degenerate eigenvalues in the ordering (|2.6[) . the 
values of cti 2 3 are non-zero according to Lemma [7J 
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Example 3 For W(x) = 1 — cosx and r\ = t/q = 0.1745, the parameters of the differential coupled-mode 
equations are approximated numerically as follows: 

ai « 0.9422, a 2 ~ 6.7813, a 3 « -4.7890. 
Let Ai,2,3(p) on p £ R 2 denote the Fourier transforms of the functions ^4i, 2,3(2/) on y £ R 2 such that 
Ai, 2 ,3(y) = 77^/ Ai,2,3(p)e ip ^, ii, 2 ,3(p) = ^ / Ai.z.aCi/Je^dj/. (4.7) 
We shall use the standard norm in L*(R 2 ) for the Fourier transforms A\p,z(p) with any q > 0: 

PIIlW) = / (l + H 2 )" /2 |i(p)Mp. (4.8) 

The differential coupled-mode equations (|4.6p are converted to the equivalent integral form 

(n - ft - a x p\ - a 2 p|) Ax(p) = aN 1 [A 1 ,A 2 ,A 3 ]{p), 

(n-P 1 -a 2 p 2 i-a 1 pl)A 2 (p) = aN 2 [A 1 ,A. 2 ,A 3 ){p), (4.9) 
(n - (3 2 - a 3 pj - a 3 pl) A 3 (p) = uN z [A u A 2 ,A 3 ]{p), 

where Ni j2 ^[Ai, A 2 , A 3 ](p) denote the Fourier transform of the cubic nonlinear terms of the differential 
coup led- mode equations. We now prove the second main result of our analysis. 

Theorem 2 Let W{x) be a bounded, piecewise- continuous, even and 2ir-periodic function on x £ R. 
Let \ < r < i. The nonlinear elliptic problem \2. 7| ) has a continuous and decaying solution $(x) for 
sufficiently small e > if there exists a solution (Bi, B 2 , B 3 ) £ L 1 (D € ,C 3 ), compactly supported on the 
disk D e = {p £ R 2 : \p\ < e r_ 2} c R 2 , which satisfies the extended coupled-mode equations in the integral 
form 

(Q-/3i- aip\ - a 2 pf) B\ (p) - oQ x (p) = e f R\ (p) , 

(n - 0! - a 2 p\ - ocxpl) B 2 {p) ~ <rQ 2 (p) = e f R 2 (p), (4.10) 
(fi - P 2 - a 3 p\ - a 3 pfj B 3 (p) - 0-Q3 (p) = e f R 3 (p) , 

where f = min(4r — 1,1 — 2r), Qi, 2 , 3 (p) denote the cubic nonlinear terms N\ t2 ^ 3 [Bi, B 2 , B 3 ](p) truncated 
on p £ D e , and Ri, 2 , 3 (p) are bounded by 



V0 < e < e : || #1,2,3 Hi* (£> e ) - Cl,2,3 ( ||-Bl|Ui(£> e ) + II-^IIl^d,) + II ^3 ||z,i( J D e ) J , (4.11) 

for some constants C\ >23 > and sufficiently small e$. Moreover, there exists an e-independent constant 
C > such that 

||$ _ _ B2 $ 2 _ B 3 $ 3 \\ c o m < Ce l ~ 2T , (4.12) 

where $1,2,3(2;) are modes defined by Lemma{^ and -Bi, 2,3(2/) are defined by the Fourier transform \j. 7[ ), 

Proof. Solutions of the nonlinear elliptic problem (|2.7p are expanded in the form 

= ^2 / $ nitn2 (ki,k 2 )u ni (xi;ki)u n2 (x 2 ;k 2 )dkidk 2 . (4.13) 

(m,n 2 )6N 2 Jj2 
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Let $ denote the union of functions ^m^i^i, ^2) f° r all ( n i> n 2) G N 2 and (£4, ^2) € T 2 and consider the 
vector space Z](N 2 ,L 1 (T 2 )) equipped with the norm 



I^II^NWr 2 )) = Yl (l + ni)%l + n 2 ) s 1^ Ij$ nun2 (k l ,k 2 )\dk 1 dk2<cc. (4.14) 

(ni,n 2 )eN 2 



T J J 



One can show similarly to the proof of Theorem Q] that the statements of Lemmas [6] and [8] extend in 
two dimensions to the vector space with the norm (|4.14p . Therefore, we convert the partial differential 
equation (|2.7p to the integral form, which is diagonal with respect to the linear terms 

[Pm(fel) +Pn 2 (k 2 ) -Uq - en]$ n (k) (4.15) 

= ~ea Yl / Mn,rn,i,j(k,l,K,X)^ m (l)^i(K)^j(X)dldKdX 
(m,i,j)eN 6 

for all n = (n\,n2) G N 2 and k = (k\, k 2 ) G T 2 , where 
Mn t m t ij(k, I, k, A) = (it ni (-; fci)^^-;^),^^-;^)^'^-; Ai))r {u n2 (-;k 2 )u i2 (-; (•;/ 2 )u J2 (-;A 2 ))r, 

for all (ni,n 2 ,mi,m 2 ,ii,i 2 , ji, j 2 ) G N 8 and (fci, fe, £1, ^2, «i 5 «2, Ai, A2) G T 8 . In view of Lemmas [6] and M 
we consider solutions of the integral equations (|4.15p in the space Z^(N 2 ,L 1 (T 2 )) for any ^ < s < 1. 

When 77 = r/o, the multiplication operator p ni (ki) + p n2 {k 2 ) — u>o vanishes at the points X' , X, and M, 
where the values of k are given by 

oi^fioVfiiHcT 2 , (4.16) 



2J ' V2 J V 2 2 . 

with n G {(1, 3); (3, 1); (2, 2)} C N 2 resp. Therefore, we apply the method of Lyapunov-Schmidt reduc- 
tions [12] and decompose the solution <E> of the integral equations (|4.15|) into 

*(fc) = E7i(fc)xDx(fe)ei,3 + U2(k)xD 2 (k)es,i + U-}(k)XD 3 (k)e 2 , 2 + (4.17) 

where {61,3,63,1,62,2} are unit vectors on N 2 , £^1,2,3 are disks of radius e r centered at the points k of 
the set (|4.16|) . see Figure [SJ Since k G T 2 , where T 2 is the first Brillouin zone, we map D on T 2 by 
periodic continuation and denote the resulting object by D n T 2 . In the representation (|4.17p . XD{k) is 
the characteristic function on k G T 2 {xD{k) = 1 if k G D n T 2 and xo(fc) = if ^ D n T 2 ), and if(k) 
is zero identically on k G Di,2 3 for the corresponding values of n: 

# ( l i3 ) =0 on D u * (3)1 ) = on L> 2 , £(2,2) = on L> 3 . 

For any of the three resonant points of A:, it follows by Lemma [7] that p' n .{kj) = and p'^ikj) 7^ for 
j = 1,2 with the n corr responding to the resonance point /c. As a result, there exists a constant C > 
such that 

min \Pm(h)\ri=rio + Pn 2 (k 2 )| ,7=770 - w | > Ce 2r . (4.18) 

(fc 1 ,fe 2 )6supp(4') 
(nj ,n 2 )eN 2 

Note that the minimal values of the multiplication operator occur for the three resonant bands with the 
indices n in the set {(1, 3); (3, 1); (2, 2)}. For all other non-resonant bands, the multiplication operator is 
continuously invertible uniformly in e for sufficiently small e. If 2r < 1, the lower bound (|4.18p is still larger 
than the perturbation terms of order e. By the Implicit Function Theorem in the space /*(N 2 , £ 1 (T 2 )) for 
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Figure 5: Decomposition of the /c-space 



any ± < s < 1, there exists a unique map * e (£>i, £> 3 ) : L^Di) x L 1 (D 2 ) x L 1 (D 3 ) ^(N 2 , L 1 (T 2 )) 
for sufficiently small e > 0. The right-hand side of the integral equation (|4.15|) is a homogeneous cubic 
polynomial of ^>{k) and Ui, 1/2,1/3 multiplied by e. Let 5 = \\Ui\\li(Di) + II^HiAflM + 
for a fixed e-independent 6q > 0. The map ^f € (Ui, U2, U3) satisfies the bound 



I^3||li(Z3 3 ) < ^0 



1*^1,^2,^3) 



U|(N 2 ,L 1 (T 2 )) 



< e x - 2r C 



\U- 



HlLi(£>i) 



+ II^2||l1(D 2 ) + ll^lki^a) ) , 



(4.19) 



for sufficiently small < e < eo, finite < 5 < 5q, and the (e, <5)-independent constant C > 0. 

We shall now map the disks -Di,2,3 to the same disk D e C M? in the stretched variable p centered at 
the origin. To do so, we apply the scaling transformation 

k — ko 



B J (p)=eU j (k), p 



=1/2 



v^Dj-nf, j = 1,2,3, 



(4.20) 



where ko is the corresponding point of the set (I4,16h , The new disk D e is now D e = {p £ R 2 : \p\ < e r ~2} 
and it covers the entire plane p 6 M 2 in the limit e — ► if 2r < 1. Note that the L 1 -norm of Uj{k) 
on k G Dj D T 2 is invariant with respect to the transformation (14.2Q[) in the new variable Bj(p) on 
p £ D e C M 2 , such that HC^-Hl 1 ^) = [l-^ilU 1 ^) f° r an y J = 1,2,3. 

When the integral equations (|4.15p are considered on the compact domains Dx, D2, D% for the three 
resonant modes and the scaling transformation (14.20P is applied to map these domains into the domain 
D e , we obtain the extended coupled- mode system (|4,10p . where remainder terms occur due to three 
different sources. The first source comes from the component * = $ e (C7i, U%, C/3) and it is estimated 
with the bound (|4.19p . The remainder terms have the order of e 1_2r and are small if 2r < 1. The second 
source comes from the expansion of p n (k) and M nm ij(k,l,K,X) in powers of e and it is estimated with 
Lemma HI The remainder terms have the order of e . The third source comes from the expansion of 
p n {k) and M n ^ m ^j(k, I, k, A) in powers oik — ko and it is estimated with Lemma[71 The remainder terms 
have the order of e 4r_1 because of the following estimate: 



HMP| 4 £j||Li(D e ) = e / H 4 B J(P) 
JD C 



dp<e 4r - 1 || J B j || L i (De) , j = 1,2,3. 



The remainder terms are small if 4r > 1. Thus, the statement is proved if \ < r < | and the largest 
remainder terms have the order e r with f = min(4r — 1, 1 — 2r). □ 
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Remark 6 If r = g, then f = r = i and both remainder terms of the extended coupled-mode equa- 
tions are of the same order in e. The fastest convergence rate of the remainder terms is e 1 / 3 because 
max re ( 1/4i i/2)(f) = 1/3. 

The differential coupled-mode system on Ai^^iu) has invariant reductions when one or two components 
of ^1,2,3 are identically zero. In particular, the reduction A3 = recovers coupled-mode equations near 
the band edge C derived in [25], while the reduction A\ = A 2 = recovers coupled-mode equations near 
the band edge B in [25]. However, these reductions do not generally persist in the extended coupled-mode 
equations near rj = r/Q since the reductions do not imply any symmetry constraints on the function <£(x). 

The coupled-mode equations (14. 9ft linearized at the zero solution describe the expansions of the function 
w = Pni(ki) + p n2 {k2) near the vertex points X' , X, and M in terms of the perturbation parameter 
e = rj — 770 and the small deviation of k = (k\, ^2) from the value ko = (kip, ^2,0) m the set (|4.16p 

w = uj + ed v iv\ v=VOtk=ko + -/^(fci.oXfci - + \p'n 2 ( k 2fi)(k 2 - h,o) 2 + 0(\k - k \ 4 ). (4.21) 

We note that the values of d v iv\ v =ri are equal for the vertex points X and X' , while the values of 
Pni(^i,o) and p'n 2 {k,2fl) have the same sign for all points X, X' and M. Moreover, sign(p'^(Xj)) = 
sign(p^. (Xj)), j = 1, 2. A simple analysis of (|4.2ip shows that the band gap between the three resonant 
Bloch band surfaces exists if 

PnSWlx > 0, p" ni {ki)\ M < 0, for ed v uj\ x > ed v u\ M 

or 

Pn 1 (ki)\x<0, p'n l (ki)\ M >Q, for ed v uj\ x < ed v uj\M- 

Examples [2] and [3] show that the first case occurs for the particular potential W(x) = 1 — cos 2 if rj > 7/0, 
such that the band gap opens for e > in the interval 02 < ^ < Pi, in a correspondence with a narrow 
band gap for 7/ > 770 in Figure E]^b) . 

5 Persistence of localized reversible solutions 

The coupled-mode system in the integral form (|4.9p is different from the extended coupled-mode system 
(|4.10j) in two aspects. First, the convolution integrals are truncated in (I4.10P on the domain D e for 
\p\ < e r ~2 : where \ < r < \. Second, the remainder terms of order e r with f = min(4r — 1, 1 — 2r) 
are present in (|4.10p . The first source is small if solutions of the extended coupled-mode system are 
considered on p £ R 2 in the space L^(R 2 ) for some q > 0. The second source is handled with the Implicit 
Function Theorem. We shall consider here a special class of decaying solutions of the coupled-mode 
system (j4.9j) on p E M 2 which lead to the corresponding solutions of the extended coupled-mode system 
(|4?TUj) on p G D t . 

Definition 1 A solution (A\, A 2 , A3) of the differential coupled-mode system (I4.6P on y G M 2 is called 
reversible if it satisfies one of the two constraints and A1.5\) for each function Ai(y), ^2(2/) an d 

Az{y). 

For notational simplicity, we write the differential coupled-mode equations (14. 6|) in the form F(A) = 0, 
where F is a nonlinear operator on (A\, A2, A3) G C 2 (M 2 , C 3 ) with the range in C°(M 2 , C 3 ). The Jacobian 
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of the operator F(A) denoted by D^F(A) is a matrix differential operator, which is diagonal with respect 
to the unbounded differential part and full with respect to the local potential part. In many cases, the 
Jacobian operator Dj^F(A) can be block-diagonalized and simplified but we can work with a general 
operator to prove the third main result of our analysis. 



Theorem 3 Let W(x) satisfy the same assumptions as in Theorem^ Let fi belong to the interior of the 
band gap of the coupled-mode system F(A) = 0. Let {A\, A2, A3) on y G R 2 be a reversible solution of the 
differential coupled-mode system F(A) = in the sense of Definition^ such that their Fourier transforms 
satisfy A G -L*(R 2 , C 3 ) for all q > 0. Assume that the Jacobian operator L>aF(A) has a three-dimensional 
kernel with the eigenvectors {d yi A,d y2 A,iA}. Then, there exists a solution (B\, B2, B3) G L 1 (D e ,C 3 ) of 
the extended coupled-mode system M-Kfy such that 

V0<e<e : \\B j - A,-|| L i (De) < C/ , Vj = 1,2,3, (5.1) 

for some e-independent constant Cj > and sufficiently small eq. 

Proof. First, we consider system (14. 10)) on the entire plane p G R 2 in the space Lg(R 2 ) for a fixed q > 
by extending the right-hand-side functions -Ri,2,3(j>) on R 2 with a compact support on D e . Thanks to 
the assumption on the existence of a solution A of F(A) = and its fast decay in L^(R 2 ) for all q > 0, 

we decompose the solution of the extended system F(B) = e r R(B) on p G R 2 into B = A + b, where b 
solves the nonlinear problem in the form 

Jb = N(b), J = D A F(A), N(t>) = e r Tt(A + b) - [f(A + b) - Jb] , (5.2) 

where J is a linearized operator, F(A+b) — Jb is quadratic in b, and R(A+b) maps an element L^(R 2 ) to 
itself. The kernel of the Jacobian operator J = Dj^F(A) is exactly three-dimensional, by the assumption, 
and it is generated by the two-dimensional group of translations in y G R 2 and a one-dimensional gauge 
invariance in arg(A). If f2 belongs to the interior of the band gap of the coupled-mode system, the 
continuous spectrum of D\F(A) is bounded away from zero and the triple zero eigenvalue is isolated 
from other eigenvalues of the discrete spectrum of DaF(A). The nonlinear elliptic problem (|1.2p with a 
real-valued symmetric potential V{x) admits the gauge invariance (fi(x) — > e ia <j)(x) for any a £ R and the 
reversibility symmetries (ll.4p and (|1.5p . The extended coupled-mode system (I4.10p obtained from the 
integral system (|4.15|) after the Lyapunov-Schmidt decomposition (|4,17p and the scaling transformation 
(|4.20p inherits all the symmetries, such that after restriction of B to functions satisfying (|1.4p or (11. 5p 
and after setting uniquely the phase, e.g. by B(0) G R, system (|5.2p is formulated in the orthogonal 
complement of the kernel of J. Therefore, the linearized operator is continuously invertible for sufficiently 
small e > under this restriction. By the Implicit Function Theorem, there exists a unique reversible 
solution in the form b = J _1 N(b) such that HbH^ip^cs) < Ce r for some C > and any q > 0. This 

result implies the desired bound (|5.ip for Bj as a solution of (|4.10p on p 6 I 2 . 

Next, we consider the error between system (I4.10p on the disk p G D e and the same system on the 
plane p G R 2 . The error comes from the terms Hb^i r D ± c s\ and ||b||xim± c 3 )' wn ere = M. 2 \D e , 



q + 2 



since A G L^(R 2 , C 2 ) for all q > and the cubic nonlinear terms Nj(b) of the coupled-mode system (|4.9 
map elements of L^(R 2 ) to elements of L^(R 2 ). Since b = J -1 N(b) and J -1 is a map from Lg(R 2 ,C 3 
to Lg +2 (R 2 ,C 3 ) for any q > 0, we obtain that 

Vb = J^Ntb) G L](R 2 ,C 3 ) ■■ ||b|Li +2(D x iC3) < llbH^^^) < C||N(b)|Ui (M2)C3) < ci 
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for some C, C > 0. Since ||b|| L i(£,± C 3) < ||b||Lj(R 2 ,c 3 ) — Ce r , the error between system (|4.1U|) on the disk 
p £ D e and the same system on the plane p S M? is of the same order as the right-hand-side terms of 
system (|4.10p . such that the desired bound (|5.ip holds on p £ D € . □ 

Corollary 2 The solution <&(x) constructed in Theorems^ and\^ is a reversible localized solution of the 
bifurcation problem {2. 7| j satisfying one of the two constraints ftl.4\ ) and 111.5]) on x 6 I 2 . 

Proof. Since the modes $12,3(20 of Lemma [2] and the solution -61,2,3(2/) °f Theorem [3] satisfy one of the 
two reversibility constraints (|1.4p and (jl.5p , the leading-order part of the representation (|4.17p produces 
a solution <&(x) which satisfies the same constraint. Since the symmetry is also preserved in the integral 
equation (|4.15p . the map $ e (C/i, U2, U3) constructed in Theorem [2] inherits the symmetry and produces 
the remainder term in the decomposition (|4.17j) satisfying the same reversibility constraint. □ 

There are several classes of localized reversible solutions of the differential coupled-mode equations 
which satisfy the assumptions of Theorem [3j These solutions are numerically approximated in Section 
[6l According to Theorem [3l all these solutions persist as localized reversible solutions of the nonlinear 
elliptic problem (jl.2p . 

Corollary 3 The solution <I>(x) constructed in Theorems^ and\^ satisfy the bound: 

||$-yli$i - A 2 <S>2 -^3^3lb°(iR 2 ) < Ce ^ f = min(4r — 1, 1 — 2r), (5.3) 

where | < r < ^. The fastest convergence rate is e 1 / 3 , which is obtained at r = ^. Consequently, the 
solution 4>{x) = y/e<&(x) of the nonlinear elliptic problem, (jl.2j) is 0{e r+l ^ 2 ) accurate to the approximation 
0( x )(x) = y^Ai^i + A2&2 + ^3^3); with the fastest convergence rate being e 5 / 6 . 

Proof. Using the triangle inequality and the bound (|4.12|) . we obtain 

3 

||$ - Ai§ 1 - A 2 $2 ~ ^3^3|lc b °(M 2 ) < Cie 1 " 2r + C 2 \\ A j - B i\\c° b {^)i 

i=l 

for some Ci, C2 > 0. By the Holder inequality we have 

\\ A j ~ ^illc^R 2 ) - \\ A j ~ -^i IIl 1 (K 2 ) = \\ A j - B j\\Ll(D e ) + ll^illLi(R 2 \-D E )' 

with the last equality due to a compact support of .61,2,3(2?) on p £ D e . The first term in the upper bound 
is 0(e r ) by the bound (|5.ip and the second term is smaller than any power of e if A 6 L*(IR 2 , C 3 ) for all 
q > 0. □ 



6 Numerical approximations of localized reversible solutions 

We approximate here several classes of localized reversible solutions of the differential coupled-mode 
system (|4.6p numerically. We use the same potential as in the Examples [EMS i.e. W(x) = 1 — cos(x) 
so that the bifurcation takes place for 770 ~ 0.1745. For selected representative cases we also verify 
the convergence rate of Corollary [3] and the conditions of Theorem [3] on the Jacobian operator, which 
guarantees persistence of these solutions in the full system (|1.2|) . We limit our attention to the following 
classes of reversible solutions: 
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(A) defocusing case a = 1: A 1 = A 2 = 0, A 3 = R(r)e m9 , r = -^=yjy\ + y§, 9 = &ig(y x + iy 2 ) 

m = ... radially symmetric positive soliton 
m = 1 ... vortex of charge one 

(B) focusing case a = —1: A3 = 



(i) A 2 = 0, Ai = fl(r)e im «, r = ^/g + g, = arg(a 2yi + *aiy 2 ) 
m = ... ellipsoidal positive soliton 

m = 1 ... ellipsoidal vortex of charge one 

(ii) A x (y x ,y 2 ) = ±A 2 (y 2 ,y x ) eR, Ai(yi,y 2 ) = Ai(-y 1) y 2 ) = A x (y x ,-y 2 ) 
. . . symmetric real coupled soliton 

(hi) A x (y x ,y 2 ) = ±iA 2 (y 2 ,yi) € iR, A 1 (y 1 ,y 2 ) = A 1 (-yi,y 2 ) = A x {y x , -y 2 ) 
. . . 7r/2-phase delay coupled soliton 

(iv) Ai(yi,y 2 ) = ±iA 2 (y 2 ,yi) G C, Ai(yi, y 2 ) = -Ai{-y x , y 2 ) = A x {y x , -y 2 ) 
. . . coupled vortex of charge one 

We will see that the localized solutions with A3 = bifurcate at a = — 1 from the upper edge Q = (3 X 
of the band gap and the localized solutions with A x = A 2 = bifurcate at a = +1 from the lower edge 
^ = /3 2 of the band gap. All solutions above satisfy the reversibility condition of Definition [TJ Provided 
the assumptions of Theorem [3] are satisfied, Corollary [2] then guarantees that these solutions correspond 
to reversible solutions of the nonlinear elliptic problem (jl.2p . Using the symmetries of ip x , ip x and (p 2 1 we 
can see that the function 4>(x) satisfies the following symmetry: 



(A) 



(B-i) 



4>(x x ,x 2 ) = (f>(-x x ,x 2 ) = (p(x x , —x 2 ) 6 R (m = 0) 
4>(x x ,x 2 ) = -<j)(-x x ,x 2 ) = 4>(x x , -x 2 ) eC (m = 1) 

(j)(x x ,x 2 ) = 4>{-x x ,x 2 ) = -(j)(x x , -x 2 ) el (m = 0) 
cf)(x x ,x 2 ) = -cj)(-x x ,x 2 ) = -</>(xi, -x 2 ) G C (m = 1) 



(B-ii) (f)(x x ,x 2 ) = ±(f)(x 2 ,x x ) = ^<f)(—x 2 , -x x ) G R 
(B-iii) 4>{x x , x 2 ) = ±i(j)(x 2 , x x ) = ^i(j>(—x 2 , —x x ) G C 
(B-iv) (f)(x x ,x 2 ) = ±i(j)(x 2 , x x ) = ±i(j>(—x 2 , —x x ) G C 

Note that (B-iii) and (B-iv) agree with the reversibility constraint (II. 5h after multiplication by e _J7r//4 . 

The above solutions do not include any three-component gap solitons with all A X ,A 2 and A3 nonzero. 
Such solutions do not bifurcate from either upper or lower edge f2 = (3 X or f2 = /3 2 , respectively, and it 
is thus hard to capture such solutions numerically. It is, nevertheless, possible that the three-component 
solitons still exist in the interior of the gap j3 2 < VL < X but we do not attempt here to locate these 
special solutions. 



6.1 One-component solutions 

One-component solutions correspond to classes (A) and (B-i). The function R(r) for class (A) satisfies 
the ODE 

R" + -R' + (O - fi 2 )R - ^R - a-f A R 3 = 0, (6.1) 
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where R(0) > 0, R'(0) = for m = and R(0) = 0, R'(0) > for m = 1. For m / the initial-value 
problem for the ODE (|6,ip is ill-posed but can be turned into a well-posed one via the transformation 
Q = r~ m R{r) leading to 

2m -I- 1 

Q" + ^—Q' + (fi - &)Q - a 74 r 2m Q 3 = 0, (6.2) 
r 

with Q(0) > 0, such that i?(r) ~ rl m l as r — > 0. Similarly, the function R(r) for class (B-i) satisfies the 
ODE 

1 TnP' 
R" + -R' + (fi - frlR - -5-i2 - (J71.R 3 = 0. (6.3) 

We solve these equations numerically via a shooting method searching for R(r) vanishing at infinity as 
r — > oo. 

Figure [6] shows the solution families in the frequency-amplitude space for both solutions with m = 
and m = 1. Class (A) families bifurcate from the lower edge O = 02 and class (B-i) families from the 
upper edge O = f3\. Examples of the gap solitons (m = 0) are in Figure [7] and of the vortices (m = 1) 
in Figure El As expected, the class (A) solutions are radially symmetric and the class (B) solutions 
ellipsoidal. 

For the case of one-component solutions, the coupled-mode equations reduce to a scalar two-dimensional 
NLS equation. Conditions of Theorem [3] on the kernel of the Jacobian operator are known to be satisfied 
for m = [28] and, therefore, do not need to be checked numerically. 

Figure [9] presents the computed e-convergence of the error term (f){x) — cj)^ l '{x) described by Corollary 
[3] for the solution classes (A) and (B-i) with m = 0. The convergence rate is seen around e 1 and e 108 , 
which is higher than the fastest convergence rate e 5 / 6 predicted by Corollary [3J The observed rate is close 
to e , which is the rate predicted by a formal multiple scales asymptotic expansion of 4>(x), for which 
(j)^(x) is the leading order term. 




Figure 6: Continuation curves of one-component solutions (A) and (B-i). The marked points P,Q, R and 
S at O, ~ 1.6, 1.22, 2.1 and 1.4 correspond to the profiles in Fig. [7]and[8l 
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Figure 7: Profiles of the one-component real gap soliton (m = 0). (a) A3 at point P in Fig. [6j (b) the 
corresponding leading-order term ^3(2/1, 2/2)^3(^1, ^2) for e = 0.1; (c) A\ at point Q in Fig. [HJ (d) the 
corresponding leading-order term A\ (2/1, 2/2)^1(^1, ^2) fo r e = 0.02. 
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Figure 8: Profiles of the one-component vortex (m = 1). (a) \A^\ at point R in Fig. [6J (b) and (c) the 
modulus and phase, respectively, of the corresponding leading-order term ^3(2/1, 2/2)^3(^1, x< i) for e = 0.1; 
(d) \A\ I at point S in Fig. [6j (e) and (f) the modulus and phase, respectively, of the corresponding leading- 
order term A\ (2/1, 2/2)^1 (^1,^2) for e = 0.1. 
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Figure 9: e-convergence of the error term for the gap soliton (m = 0) corresponding to (A) at ~ 1.22 
and (B-i) atfi« 1.25. 

6.2 Two-component solutions 

The following algorithm is used to compute two-component solutions of classes (B-ii) , (B-iii) and (B-iv) . 
First, a\ and «2 are replaced by a = ^(01+02) so that the first two equations of the coupled-mode system 
(I4.6P admit a reduction A± = A2, where A\ = R{r)e imd solves a one-component problem. The function 
R(r) is then computed as in Section 16.11 Next, homotopy in the coefficients a± and 0:2 is employed to 
gradually deform the computed solution with a\ = a.2 = a into the one with the original values of a\ and 
a2. In this homotopy continuation the full coupled- mode system needs to be solved numerically and this 
is done via Newton's method on a central finite-difference discretization with zero Dirichlet boundary 
conditions. 

Figure [10] shows the solution families of classes (B-ii) , (B-iii) and (B-iv) as curves in the frequency- 
amplitude space. Profiles of examples of these solutions corresponding to the marked points in Figure [TOI 
appear in Figures [Til E2] and [13] Note that the coupled vortex in Figure [13] has been obtained via the 
above described homotopy continuation from a vortex of charge one with a\ = «2- As the phase plots of 
A\ and A2 show, the resulting coupled vortex is also of charge one. 

In order to verify the persistence conditions of Theorem [3] for the two-component solutions, the kernel 
of the Jacobian operator J has to be studied numerically. When rewritten in the real variables (for real 
and imaginary part) the first two equations of the coupled-mode system (|4.6p give rise to a Jacobian 
operator J. For the cases (B-ii) and (B-iii), the Jacobian operator can be block-diagonalized to the form 



such that dim(Ker(J)) = dim(Ker( J + )) + dim(Ker( J_)). The elements of Ker(J) can, moreover, be easily 
reconstructed from the elements of Ker(J + ) and Ker(J_). This diagonalization reduces the expense of 
the eigenvalue solver. 

Figure [141 shows the four smallest (in modulus) eigenvalues Xj ± , . . . , Aj 4 of the 4th order central finite- 
difference approximation of J for examples of solutions (B-ii) and (B-iii) as functions of D, where D is 
the size of the computational box [— D, D] x [— D, D] C M? with the step sizes dyi = dy2 ~ 0.14 and with 
zero Dirichlet boundary conditions. Clearly, the three smallest eigenvalues converge to zero as D grows 
while the fourth eigenvalue converges to a nonzero value. Therefore, the assumptions of Theorem [3] are 
verified for these solutions. We have also checked that the corresponding eigenvectors of the three zero 
eigenvalues are approximations of d Vl A,d y2 A and iA. 
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For the case (B-iv) the Jacobian operator J has to be treated in full because it cannot be block- 
diagonalized for complex- valued solutions. Figure [15] shows the four smallest eigenvalues of the discretized 
Jacobian operator for an example of the solution (B-iv). The fourth-order central finite-difference dis- 
cretization was used once again with the step sizes dy\ = dy 2 ~ 0.12. The fourth eigenvalue A j 4 converges 
approximately to 0.005 as D grows, while the second and third eigenvalues Aj 2 3 converge to 0.0002. We 
have checked that the limit of Xj 23 decreases as dy\ = dy 2 is decreased while the limit of Aj 4 remains 
practically unchanged. 

Finally, Figure [16] verifies the e-convergence of the error term (j)(x) — (fP^ (x) described by Corollary 
[3] for the gap soliton of class (B-ii) at ~ 0.944. The observed convergence rate is e 95 , which is, 
once again, close to e 1 expected from the formal asymptotic expansion. The e-convergence has not been 
checked for any complex-valued solutions of classes (B-iii) and (B-iv) due to four times larger memory 
requirements compared to the real case. The memory use grows rapidly as e decreases since domains of 
the size (2D/y / e) x (2D/y/e) need to be discretized to compute the solution 4>(x) of the elliptic problem 
(jl.2p and each period [0, 2tt] x [0, 2n] of the potential function V(x) requires about 100 grid points. 




Figure 10: Continuation curves of two-component solutions (B-ii), (B-iii) and (B-iv). The marked points 
T and U at Q ~ 2.04 correspond to the profiles in Figs. [11] and [T2"j while the point V at f2 ~ 1.4 
corresponds to the profiles in Fig. [13] 
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Figure 11: Two-component real gap soliton (B-ii) at point T in Fig. [10] (a) A\\ (b) A 2 ; (c) the leading- 
order term A 1 (y 1 ,y 2 )®i(x 1 , x 2 ) + A 2 (yi, 2/2)^2(^1, x 2 ) for e = 0.1. 
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Figure 12: Two-component 7r/2-phase delay gap soliton (B-iii) at point U in Fig. [TU1 (a) A\\ (b) 
—%Ai\ (c) and (d) modulus and phase, respectively, of the leading-order term Ax(yx, U2)^i(xi, x 2 ) + 
M{yi, 2/2)^2 (xi, x 2 ) for e = 0.1. 
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Figure 13: Two-component vortex (B-iv) at point V in Fig. [TUJ (a) and (b) modulus and phase of 
Ai; (c) and (d) modulus and phase of A2; (e) and (f) modulus and phase of the leading-order term 
Ai(yi, 2/2)^1(^1,^2) + A 2 (y 1 , 2/2)^2(^1, x 2 ) for e = 0.1. 
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Figure 14: The four smallest eigenvalues of the Jacobian operator J for (a) the two-component real gap 
soliton (B-ii) and (b) the two-component 7r/2-phase delay gap soliton (B-iii) for Q ~ 1.19. 
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Figure 15: The four smallest eigenvalues of the Jacobian operator J for the two-component vortex (B-iv) 
at n 1.2. 




Figure 16: e-convergence of the error term for the gap soliton (B-ii) with f2 ~ 0.94. 

7 The time-dependent case 



We show here that dynamics of non-stationary localized solutions of the Gross-Pitaevskii equation (jl.ip 
are close to the dynamics of the time-dependent coup led- mode equations for finite time intervals. These 
results are similar to the ones in [23] and they give a rigorous basis for formal asymptotic results in 
[6j[25]. According to a formal asymptotic multiple scales expansion, solutions of the Gross-Pitaevskii 
equation (jl.ip in the form 



E ans (x,t) = ^E 1 (x,t) + 0(e), 



where 



Ei = [Ax (\fex, et)ipi {xx)(p 2 {x 2 ) + A 2 (y/ex, et)cp 2 (xx)i^x(x 2 ) + A 3 (y/ex, et)<px (xx)<Px(x 2 )] e 
are approximated by solutions of the time-dependent coupled-mode equations 



(id T - (3 1 )A 1 

{id T - Pi)A 2 + (a 2 d 2 yi + a 1 9i 2 )A 2 
(id T - f3 2 )A 3 + a 3 (dl+dl 



a 2 dl)A x 

2 S 



A, = 



Jl \A 1 \ 2 A 1+l2 (2\A 2 \ 2 A 1 +AlA 1 ) 
ll \A 2 \ 2 A 2+l2 (2\A 1 \ 2 A 2 + A 2 A 2 ) 



. 7 3(2|^ 3 | 2 ^i+^Pi) 
■ l3 (2\A 3 \ 2 A 2 +AlA 2 ) 



14 \A 3 \ 2 A 3 + 2 l3 (\A 1 \ 2 + \A 2 \ 2 )A 3 + 13 {A\ + A 2 )A 3 



where y = ^fex and T = et. Under the non-resonance condition 



inf 

(ni ,n 2 )eN 2 
J1,J2.J36{±5,±3,±1} 



h + h 



+ Pn 2 



31 + 33 



\jl +32 + J3 1 ^0 



>0, 



(7.1) 



(7.2) 
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where IJ1+J2+J3I < 5 and the infinimum does not include direct resonances that give the nonlinear terms 
in system (|7.ip . the following theorem states the approximation result on the time-dependent solutions. 

Theorem 4 Let s > l,s* > max{3, s} ; and assume that the non-resonance condition (|7.2p holds. 
Then for all C\ and Tq > i/iere exist eo > and C2 > smc/i i/iaf /or a// solutions A\,A2, A3 E 
C([0,Tq], H s * (M?)) of the time- dependent coupled-mode system with 

sup \\A(',T)\\ Ha ,, R 2) < Ci, 

Tg[0,T ] 

there exist solutions E E C([0, To/e], H S (M?)) of the Gross-Pitaevskii equation (jl.ip tmt/i 

sup ||£(.,t) -Kns(-,i)||^(R2) < C 2 e 3 / 2 . 

*6[0,T /e] 

/or a// e E (0, eo). 

Proof. The proof is very similar to the one given for a scalar nonlinear Schrodinger equation in the 
one-dimensional case in [7J. The non-resonance condition (|7.2p differs from Eq.(4)] by including the 
terms with ji,j2)j3 = ±5. This is necessary to make the residual formally of order 0(e 3 ). Due to the 
fact that we are working here on d 6 R 2 , we lose e -1 / 2 due to the scaling of the L 2 -norm in contrast to 
[7J, where we would only lose e" 1 / 4 on x E M. with the same scaling. Assuming that the non-resonance 
condition (|7.2j) is satisfied, the 0(e) terms in .E ans can be chosen in such a way that all terms up to order 
0(e 5 / 2 ) in the residual 

Res(£ ans ) = -%Et - V 2 E + V(x)E + a\E\ 2 E 

can be eliminated. By [7J Lemma 3.3] Bloch transform is an isomorphism between H S (R 2 ) and l 2 s (N 2 , L 2 (T 2 )). 
Hence, similarly to Lemma 3.4] the nonlinearity E t— > \E\ 2 E maps / 2 (N 2 , L 2 (T 2 )) in Bloch space, or 
iF(R 2 ) in physical space, into itself for s > 1 due to Sobolev's embedding theorem in two dimensions. 
Moreover, L = —V 2 + V(x) generates a uniformly bounded semigroup in / 2 (N 2 , -L 2 (T 2 )), or respectively 
in H S (M 2 ), according to [TJ page 927]. Hence, the error function R = e~ 3 / 2 (E — E ans ) satisfies an equation 
of the form 

d t R = LR + 0(e). 

Using the variation of constant formula and Gronwall's inequality gives the required 0(1) bound for R. 
For more details see [7j Section 4.2]. □ 

8 Various generalizations 

We have given a detailed justification of the three coupled-mode equations which describe gap solitons 
bifurcating in the first band gap of a two-dimensional separable periodic potential. These equations are 
generic (structurally stable) for the considered bifurcation but they can be modified for other relevant 
bifurcation problems. We review here several examples of other bifurcations, where the justification 
analysis is expected to be applicable in a similar manner. 

Bifurcations from band edges. If the band edge separates the existing band gap of finite length and a 
single or double-degenerate spectral band, the coupled-mode equations persist without any modifications. 
A formal derivation is given in the recent paper [25], where a scalar NLS equation is derived for a single 
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band and two coupled NLS equations are derived for a double band. Note that the coupled NLS equations 
for a double band may have terms destroying the reduction A\ = or A 2 = in the truncated system of 
equations. This happens ,for instance, at the band edge E in the second band gap of a separable symmetric 
potential [25J. Since the analysis of separable potentials relies on the construction of one-dimensional 
Bloch modes and one-dimensional spectral bands, the same analysis is valid for bifurcations from band 
edges in one-dimensional problems. The formal derivation of the NLS equation for a band edge of a single 
spectral band in the one-dimensional GP equation was reported in |22j . 

Bifurcations in the higher- order band gaps. With larger values of the parameter 77 in the separable 
potential (jl.3p . more band gaps open, although the number of band gaps is always finite for finite 77, see 
Theorem 6.10.5 in [10]. Since the bifurcation of new band gaps occurs due to the resonance of finitely 
many Bloch modes in the separable potential, a system of finitely-many coupled-mode equations can be 
derived similarly in the higher-order band gaps. For instance, the bifurcation of the second band gap in 
the symmetric separable potential may occur for a resonance of either three modes when A2 + ^1 = 2/i 2 < 
Ai + A3 or four modes when A2 + ^1 = Ai + A3 > 2^2 or five modes when A2 + fJ-i = Ai + A3 = 2/1/2- 

Bifurcations in asymmetric separable potentials. If the separable potential is asymmetric, such that 

V{x)=r ] (W 1 {x 1 ) + W 2 (x 2 )), 

where W\ and W 2 are two different functions on x € R, the degeneracy of spectral bands is broken and 
bifurcations of the band gaps occur with a smaller number of resonant Bloch modes. Our analysis remains 
valid but the spectrum of two different one-dimensional Sturm-Liouville problems must be incorporated 
in the Fourier-Bloch decomposition. In this case, two coupled-mode equations occur generally for the 
bifurcation of the first band gap in the anisotropic separable potential. 

Bifurcations in finite-gap potentials. Finite-gap potentials lead to the degeneracy of eigenvalues of the 
one-dimensional Sturm-Liouville problem and non-zero derivatives of the functions p n (k) at the extremal 
points k = and k = ±i. Adding a generic potential breaks the symmetry of the finite-gap potential and 
leads to the bifurcation of narrow band gaps in the space of one dimension. The corresponding coupled- 
mode equations must have first-order than second-order derivative operators. A formal derivation of such 
coupled-mode equations from the Bloch mode decomposition was considered in [9]. Our analysis provides 
a rigorous justification of these coupled- mode equations. 

Bifurcation in super-lattices. Let the potential V{x) be represented in the form 

V(x) = 7 ?0 (W( X1 ) + W(x 2 )) + e (W(x!) + W(x 2 j) , 

where W{x) and 770 are the same as in the separable potential (|1.3p and W[x) is a 47T-periodic potential. 
The perturbation term W(x) couples resonant Bloch modes of the potential W{x) by linear terms and 
destroys reductions A\ = 0, A 2 = and A3 = in the differential coupled-mode system. The system is 
still formulated in the form of the coupled NLS equations with various (linear and nonlinear) coupling 
terms. 

Bifurcations in three-dimensional periodic problems. The same analysis holds in three-dimensional 
separable potentials, since the spectral bands are still enumerated by a countable number of spectral 
bands of the one-dimensional potentials. Bifurcations of new band gaps occur again due to a resonance 
of finitely many Bloch modes of the one-dimensional potentials. 

Bifurcations in non-separable periodic potentials. At the present time, the generalization of the analysis 
for non-separable periodic potentials meets a technical obstacle that the bound (|3.7f) obtained in [7] is 
needed to be extended to problems with two-dimensional periodic potentials. 
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